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基于 旋转 基 技 术 的 多 目标 粒子 群 优化 算法 
戴 永 彬 ， 纪 川 川 ， 康 伟 伟 


(辽宁 工业 大 学 软件 学 院 , 辽宁 锦州 121001) 


摘 要 : 为 了 解决 多 目标 优化 求解 的 问题 ， 提 出 一 种 基于 旋转 基 技 术 的 多 目标 粒子 群 优化 算法 (rtMOPSO)。 首 先 ,改进 

了 旋转 基 可 视 化 技术 并 将 Pareto IS RH AMEN RHR HBP HL, RA AWA ARM AeA; 其 次 ， 为 平 
衡 归 档 集 的 收 么 性 和 多 样 性 ， 提 出 了 角度 支配 和 角度 支配 力度 两 种 新 的 概念 并 设计 归档 集 新 的 排序 方法 ; 最 后 ， 在 融 

合 了 旋转 基 角 度 和 距离 的 概念 的 基础 上 ， 提 出 了 一 种 改进 的 全 局 引导 粒子 的 选择 策略 。 改 进 算 法 采用 两 个 类 别 的 测试 

函数 ， 与 $ 种 多 目标 优化 算法 进行 了 对 比 实 验 。 实 验 结果 表 明 ， 改 进 算法 在 收敛 性 和 多 样 性 方面 优势 明显 。 
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Multi-objective particle swarm algorithm based on rotation basis 


Dai Yongbin, Ji Chuanchuan, Kang Weiwei 
(College of Software Liaoning University of Technology, Jinzhou Liaoning 121001, China ) 


Abstract: This paper proposed a multi-objective particle swarm algorithm based on rotation basis (TtMOPSO) in order to solve 
problems about multi-objective optimization. Firstly, Pareto front maps to the improved sectorial plane of rotation basis so that 
evolutionary status of the population were detected by the entropy and its difference entropy of the population. Then, angle 
dominance and strength of angle dominance were proposed to design an archive maintaining strategy which can balance diversity 
and convergence. Finally, selecting the global best solution was proposed based on rotation angle and distance of the rotation 
basis. The improved algorithm is compared with 5 multi-objective optimization algorithms on two kinds of test suites. The 
Simulation results show that improved algorithm in this paper has big advantages over the other competitors in terms of diversity 


and convergence. 
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可 以 采取 措施 避免 粒子 种 群 早熟 ， 但 是 一 定 会 影响 粒子 种 群 后 
期 的 多 样 性 。 因 此 如 何平 衡 局 部 收敛 性 和 种 群 多 样 性 急需 有 效 
近年 来 ,多 目标 优化 问题 受到 了 学 术 界 和 工程 界 广泛 关注 ， 的 解决 方案 。 为 此 ， 进 一 步 完善 归档 集 的 管理 与 维护 、 全 局 引 
取得 了 大 量 的 研究 成 果 。 然 而 当 多 目标 优化 问题 的 维 数 大 于 3 导 粒 子 的 选取 以 及 种 群 进化 状态 的 监测 等 方面 的 算法 成 为 提高 
时 ， 如 何 准 确 呈 现 Pareto 前 沿 的 分 布 状 态 从 而 实施 相应 的 控制 。” 粒子 群 收敛 性 和 多 样 性 的 关键 所 在 。 
策略 获取 最 优 Pareto 解 集 ， 是 一 项 十 分 困难 的 工作 。 可 视 化 技 前 ， 在 归档 集 的 管理 方面 ， 国 内 外 学 者 提出 了 多 种 维护 
术 是 一 种 解决 多 目标 优化 问题 的 重要 手段 种 。 可 视 化 技术 主要 ”策略 进 ， 主 要 分 为 两 类 : 第 一 类 方法 是 根据 种 群 密度 的 大 小 剪 
采用 降 维 算法 映射 Pareto 前 端 ， 分 为 无 损 降 维和 有 损 降 维 。 无 ，” 裁 归档 集 ， 例 如 ， 自 适应 网 格 、 拥 挤 距离 、6 支配 外 等 方法 。 
损 降 维 方法 实现 降 维 不 会 损坏 数据 信息 , 主要 包括 平行 坐标 系 、 ”第 二 类 方法 是 依据 归档 集 粒子 和 参考 点 或 参考 区 的 距离 整理 归 
旋转 基 坐 标 等 方法 9; 有 损 降 维 首先 获取 高 维 数据 的 特征 和 关 ” 档 集 ， 例 如， 参考 区 域 、 参 考点 等 偏好 算法 9% 由。 在 全 局 引导 
系 ， 将 其 映射 到 低 维 空间 ， 主 要 包括 主 成 分 分 析 、 高 斯 随机 过 ” ”粒子 的 选取 方面 ， 现 已 提出 了 大 量 的 性 能 优异 的 算法 和 改进 算 
程 模型 、 聚 类 等 659。 文献 [和 ] 利 用 旋转 基 可 视 化 技术 反映 Pareto 法 , 如 sigma 以 及 改进 算法 、 灰 色 关 联 、 信 息 烂 法 、 差 分 算法 、 
解 优 劣 情况 ， 方 便 决 策 者 进行 分 析 和 决策 ， 但 是 不 能 监控 粒子 。” 动态 加 权 、 轮 盘 赌 I> 咏 等 方法 。 通 过 这 些 算法 可 以 获得 全 局 最 
种 群 的 进化 状态 ， 有 效 解决 多 目标 优化 问题 。 优 粒子 ， 从 而 保证 了 粒子 种 群 的 正确 飞行 方向 ， 提 高 了 种 群 的 
粒子 群 可 以 快速 收敛， 但 容易 早熟 ， 陷 入 局 部 最 优 。 虽 然 。 ”收敛 性 和 分 布 性 。 另 外 ， 还 有 多 种 根据 种 群 进化 状态 调节 收敛 
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性 和 多 样 性 的 平衡 算法 ， 如 动态 调整 惯性 权重 和 学 习 因 子 、 种 ”数据 挖掘 、 网 络 、 社 交 、 文 本 等 可 视 化 方面 。 旋 转 基 技术 是 重 
群 规模 自 适 应 算法 、 种 群 的 交叉 与 变异 算法 、 多 层次 信息 交互 ”要 可 视 化 技术 之 一 ,可 以 将 多 维 数据 直接 呈现 给 用 户 或 决策 者 ， 
算法 等 ， 但 是 这 些 算法 对 种 群 进化 状态 还 未 能 准确 监测 ， 因 此 甫 助 决策 。 但 该 算法 不 适合 直接 用 于 种 群 进化 状态 的 监控 ， 文 
也 就 不 能 清楚 地 辨别 进化 状态 ， 从 而 有 效 地 执行 进化 策略 。 为 ” ” 献 [和 也 没有 将 其 用 于 多 目标 优化 问题 的 求解 。 为 此 ， 本 文 对 旋 
了 监测 种 群 进化 状态 ， 文 献 [8] 提 出 了 状态 烂 以 及 差 灶 的 概念 ， ” 转 基 进 行 了 改进 并 应 用 于 多 目标 粒子 群 优化 算法 。 


对 归档 集 状 态 糯 的 变化 情况 进行 分 析 和 辨识 ， 从 而 确定 粒子 群 2.1 旋转 基 技 术 的 改进 

进化 状态 的 边界 阔 值 。 但 文献 [8] 使 用 密度 排序 方法 维护 归档 集 ， 旋转 基 可 视 化 算法 的 基本 原理 是 利用 旋转 基 向 量 的 旋转 f 

这 种 单一 排序 方法 不 能 同时 满足 归档 集 收敛 性 和 多 样 性 的 维护 ERR Pareto 解 在 各 维 目标 上 的 取 值 情况 ， 再 将 各 目标 上 的 基 

需求 。 a), 在 90 扇形 平面 上 形成 一 条 折线 。 绘制 点 的 旋转 角度 
综合 以 上 分 析 ， 本 文 将 平行 坐标 和 文献 [和 1] 提出 的 旋转 基 技 ” 表示 该 解 的 性 能 ， 终 点 向 量 到 原点 的 长 度 或 距离 表示 该 解 性 能 

术 相 结合 ， 改 进 了 旋转 基 映 射 策略 从 而 实现 高 维 数据 到 二 维 空 。” 的 波动 情况 。 标 准 化 后 的 旋转 角度 公式 如 下 所 示 : 

间 的 映射 ， 使 决策 者 可 以 对 Pareto 解 集 的 进行 准确 地 分 析 和 决 Pe ied 人 (3) 

策 。 其 次 ， 在 归档 集 的 维护 和 管理 方面 提出 了 旋转 基 角 度 支 配 fR 2 

和 角度 支配 力度 的 新 概念 ， 设 计 了 归档 集 双 排序 算法 。 另 外 ， 其 中 : i 为 目标 的 数量 ， 二 1,2，…M; 大 表示 种 群 的 粒子 的 

本 文 将 旋转 基 角度 和 距离 作为 全 局 最 优 粒子 的 选择 指标 ， 因 此 ” 序号， 本 1,2，… Ks 户 、 记 为 适应 度 函 数 的 最 大 值 和 最 小 值 。 

可 以 获取 综合 性 能 最 强 的 粒子 作为 全 局 最 优 粒 子 。 最 后 ， 对 文 旋转 基 向 量 的 连接 过 程 如 下 : 首先 ， 设 定 一 组 固定 长 度 的 


献 [8] 提 出 的 基于 Pareto 粒 的 状态 监测 方法 进行 了 改进 , 通过 监 PEE (wi, 计 1,2,… M}, 向 量 的 数量 等 于 目标 的 维 数 。 其 
测 种 群 状态 选择 不 同 的 归档 集 排序 方法 ， 实 现 了 归档 集 的 综合 次， 以 原点 为 起 点 ， 按 照 水 平方 向 ， 逆 时 针 旋 转角 度 % ， 得 到 
动态 管理 ， 提 高 了 种 群 的 收敛 性 和 分 布 性 。 Pareto 解 的 第 一 个 目标 的 映射 坐标 , 然后 再 以 该 点 为 起 始 坐标 ， 


日 ` X) J 同 单方 式 绘 Hl) 下 不 WN 的 映射 ZEER o 依次 类 推 即 J 完成 
1 多 目 标 优 {4 > 问题 和 米 J f 群 基本 原理 
第 M 个 目标 映 稍 Al fr Fk 绘制 。 


= 
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1.1 多 目标 优化 问题 综 上 分 析 可 知 ， 旋 转 基 可 视 化 技术 可 以 表示 Pareto 解 的 性 
多 目标 问题 可 描述 如 下 : 能 的 优 劣 程度 ， 但 是 没有 采用 类 似 密度 等 指标 展示 全 部 粒子 在 


eee ee Pareto 前 沿 分 布 状况 四。 为 此 , 本 文 将 旋转 基 和 平行 坐标 结合 起 
r 来 ， 改 进 旋转 基 技术 。 有 具体 如 下 : 根据 种 群 粒子 的 数量 将 旋转 
Eee FE O ” 基 90 度 遍 形 区 域 划 成 天 个 扇 区 ， 每 个 扇 区 的 圆心 角 的 角度 为 
到 /2K。 然 后 采用 平行 坐标 的 思想 ， 将 目标 维 数 取 M 条 垂直 凶 
线 ， 每 条 轴线 上 的 数值 范围 为 目标 适应 度 函 数 的 最 小 值 到 最 大 
值 并 且 均 匀 分 布 。 为 了 将 旋转 基 融 入 平行 坐标 ， 本 文 将 平行 坐 
标的 M 条 垂直 轴线 转换 为 旋转 基 90 度 扇 区 内 的 M 条 同心 圆 圆 
弧 , 这 样 就 形成 KxM 个 扇 块 。 经 过 平行 坐标 映射 的 旋转 基 平 本 
称 为 旋转 基 扇 区 平面 ， 被 同心 圆 圆 弧 所 分 割 的 扇 区 平面 称 为 旋 
HEIER. W, $ K=3，M=4， 这 样 形成 共计 12 个 扇 块 ， 有 具体 
如 图 1 所 示 。 


其 中 : XOW R 7ST RRA g(x) AA 万 0 分 别 是 不 等 式 约 
束 和 等 式 约 束 。 

多 目标 优化 问题 一 般 采 用 Pareto 解 集 作 为 寻 优 的 结果 。 
Pareto 解 集 内 一 般 有 若干 非 劣 解 ， 可 以 根据 决策 需要 选取 合适 
的 非 劣 解 作为 最 优 解 ， 在 Pareto 前 沿 均匀 分 布 并 无 限 靠近 真实 
值 ， 这 是 Pareto 解 集 的 最 理想 的 状态 。 

1.2 ”粒子 群 算法 基本 原理 

粒子 群 是 一 种 基于 群体 仿生 的 进化 算法 。x 代表 粒子 的 位 
B, v 代表 粒子 的 速度 ， 两 者 的 关系 可 以 描述 为 


rH t t t t t 
v” = wv; +anCp; — x) +c, lg; — x) 


到 


SS 


(2) 


x =x tit 
其 中 : v S= OV Vip) X =O Xin Xp) Ah, w 为 惯性 
权 值 ， 个 体 最 优 值 和 全 局 最 优 值 分 别 是 
Pi=(pins Post Po) 8 =(81,82,…,8p)，c1、C2 是 学 习 因 子 ， 
ri， 允 是 介 于 [0,1] 的 随机 数 ，t 表示 粒子 群 迭 代 更 新 次 数 ，D 


代表 粒子 的 维 数 。 A1 旋转 基 肩 区 平面 
这 样 ， 高 维 的 Pareto 前 沿 被 映射 成 旋转 基 扇 区 平面 内 ， 种 
2 ”旋转 基 技术 及 进化 状态 监测 的 改进 群 粒子 也 同时 被 映射 到 旋转 基 扁 块 内 。 通 过 分 析 遍 区 块 内 的 粒 


子 的 数量 可 以 获知 Pareto 前 沿 粒 子 的 分 布 状态 。 为 了 计算 每 个 
粒子 的 目标 分 量 所 在 的 扇 块 位 置 ， 可 用 如 下 公式 表示 。 如 第 无 


多 维 数据 可 视 化 属于 信息 可 视 化 的 研究 领域 ， 主 引 
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录用 稿 


个 粒子 第 i 个 目标 分 量 映射 到 旋转 基 扇 区 的 扇 块 编号 如 为 


Ea K-06 
nee ° 


其 中 : [x] dai LM, APN, KIB, mR OF WO 
时 , AAS A 1. Pl, BEL A eit ial A 4 个 目标 ， 粒 
子 总 数 为 3， 按照 公 式 (4), 可 计算 出 每 个 粒子 编号 。 例 如 ,第 
1 个 粒子 的 第 2 个 目标 分 量 编号 为 3， 即 hi 的 五 值 为 3, 可 写 
成 有 ;=3。 具 体 分 布 如 图 2 所 示 。 


图 2 ”旋转 基 扇 区 平面 中 Pareto 前 沿 
2.2 IK RHMKA HA 
2.2.1 RASA 

本 文采 用 了 Pareto Wi FZ OUT WY Pareto 前 沿 的 进化 状态 。 
Pareto 烂 和 差 和 能 够 准确 描述 旋转 基 扇 区 平面 上 种 群 粒子 分 布 
情况 。 

1) REAR 


àT 


E(t) = 


£ W Sector” (t), Sector” (t) 

-之 之 KM °° KM (5) 
其 中 : Sector! CD) 表示 落 在 第 HATH i IARA KER 

量 的 数量 。 
2) FRAR 


1 


AE = E(t) -Elt —1) (6) 

2.2.2 进化 状态 阅 值 的 确定 
粒子 群 种 群 在 进化 过 程 中 会 经 历 收敛 状态 、 多 样 性 状态 和 
停滞 状态 等 三 个 状态 。 其 中 ， 多 样 性 状态 为 种 群 收 依 和 停滞 的 
中 间 状 态 ， 因 此 只 需要 确定 中 间 状 态 差 粒 的 最 大 边界 冰 值 a 和 
最 小 冰 值 6， 就 可 以 确定 进化 状态 。 为 此 ， 设 定 种 群 分 布 的 一 
种 极端 情况 : Pareto 最 优 解 集 有 3 个 粒子 且 有 3 个 目标 分 量 
粒子 k 的 所 在 的 每 个 户 块 内 都 有 两 个 分 量 ， 有 具体 如 图 3 所 示 。 
当 Pareto 前 沿 迁 代 更 新 后 ， 出 现 的 第 一 种 极端 情况 是 ， 旧 粒子 
k WEIT k 替换 并 全 部 落 入 了 3 个 空闲 的 扇 块 ， 有 具体 如 
图 4 ; 另 一 种 情况 是 更 新 粒子 局 只 落 入 了 1 个 空闲 的 扇 块 。 归 
MARA RAT MN, HERRER, HRA EI 
RLF FEMA AR, KER SEP EH ET BI iw AER 
AY ete BN Ay i VSS Se ee BARUR (7) (8) 所 示 。 

«(San e Lio bea) 

3 (7) 


m i 1 40 
(= 而 log om 2 oM log ul ae 


oie 
B= UE log2 (8) 


该 根据 种 群 进化 状态 
了 保持 收敛 速度 ,保留 收 敛 性 强 的 粒子 ， 当 
态 时 ， 为 维护 归档 集 的 多 样 性 ， 删 除 分 布 密 
为 此 ， 本 文采 用 双 排 序 的 维护 策略 ， 在 不 同 的 进化 阶段 选择 不 
同 排序 策略 ， 平 衡 选 择 压力 和 分 布 性 。 本 文 提 HH 
标 和 分 布 程序 指标 作为 双 夺 
3.1.1 收 你 程度 指标 

为 了 度量 种 群 的 收敛 程度 ， 本 文 提 出 了 角度 支配 和 角度 支 
配 力度 的 概念 并 将 角度 支配 力度 作为 收敛 程度 的 指标 。 


BOO, Kaxi 


其 中 ， 为 O SHAMS A 4 


归档 集 的 维护 与 全 局 引导 粒子 选取 


归档 集 的 维护 与 管理 


粒子 群 的 全 部 非 支配 解 一 般 存放 在 ) 
归档 集 维护 策略 是 实现 归档 集 


多 数 归 档 集 采 用 单一 的 管理 策 
单 性 指标 。 例 如 ， 归 档 集 经 常 ? 
指标 加， 这 种 策略 侧重 于 


AL, HAE Pa È OOOG-1) 时 ， 种 群 处 于 收敛 状态 ; 
当 4=K 且 p<|AE |<a 时 处 于 多 样 化 状态 ; 4 A=K HAE [<L 时 ， 
种 群 处 于 停滞 状态 。 


图 3 更 新 前 的 情况 


图 4 更 新 后 的 情况 


选择 压力 的 控制 ; 


种 则 侧重 于 分 布 性 的 控制 。 


实际 上 ， 归 档 集 的 


自 适应 切换 ， 当 种 群 处 于 收敛 状态 时 ， 为 


度 比 较 高 的 粒子 。 


1) 角度 支配 概念 


归档 集 的 任意 两 个 解 为 x F y, P 
Exi y SAMS 
Vi=1,2,--M,0 <6 AFj=1,2,--M,0! <0 


F 序 的 依据 ， 
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2) 角度 支配 力度 概念 
归档 集 的 任意 一 个 解 xe4， 其 支配 力度 Si 是 角度 支配 其 
他 解 的 数量 。 


| 


Si={u|x EA, *~<u } (10) 
角度 支配 力度 表明 了 归档 集中 粒子 向 Pareto 前 沿 收敛 的 程 
度 ， 角 度 支 配 力度 越 大 ， 说 明 该 解 越 靠近 真实 的 Pareto 前 沿 。 
3.1.2 分 布 程度 指标 


eee 
F= JH aa) (11) 
其 中 : K X Pareto 解 集中 解 的 总 数 ，; 
d, =min(》 Hi (x) Hil) ,， Ciz7) 表 示 某 个 非 劣 解 的 目标 向 


到 最 优 解 目 


ty 


BARRER: JÆ di 的 平均 值 。 


ASS 


起 。 


保证 


排序 策略 .总 所 周知 ,密度 指标 
因此 ， 该 策 


和 E 护 和 全 局 最 优 粒子 选取 进行 了 
KRF. He CH 


改进 ， 使 本 文 算法 的 性 能 有 较 


省 没有 将 粒子 排 


[8] 对 归档 集 的 允 


本 文 算 法 可 根据 种 群 的 进化 状态 ， 采 | 
了 归档 集 粒 子 的 多 样 性 和 收敛 性 。 


EPA} 


密度 指标 进行 单 
主要 反应 归档 集 粒子 的 多 样 


ti 


序 情况 和 种 群 的 收敛 性 能 结合 在 一 


表 1 格 占 优 强度 与 旋转 基 角 度 


j 双 排序 策略 ， 同 时 


分 布 程度 指标 Fi 反应 了 非 劣 解 集 对 应 前 端的 均匀 程度 。 划 
值 越 小 ， 分 布 性 能 越 好 。 
3.1.3 旋转 基 折 线 距 离 
归档 集 粒 子 排序 时 ， 如 果 粒 子 S 和 F 的 指标 都 相同 ， 则 
用 旋转 基 上 映射 折线 距离 作为 辅助 判断 指标 。 旋 转 基 折 线 距 离 
公式 ， 具 体 如 下 : 


a 


a (12) 
旋转 基 折 线 距 离 越 大 说 明 该 粒子 在 各 个 目标 上 性 能 波动 越 


小 ， 性 能 越 均衡 。 
3.1.4 归档 集 的 排序 与 维护 

根据 进化 状态 ， 如 果 种 群 处 于 收敛 阶段 ， 先 按照 St、Fi 和 
Lt 次序， 对 关键 字 5; 按 升序 而 Fi 和 Lx 按 降序 排序 ， 如 果 种 群 
处 于 多 样 化 阶段 ， 则 归档 集 粒子 可 按照 Fi Si 和 Lx 次序， 对 对 
关键 字 8 按 升序 而 F; 和 Li 按 降序 排序 。 最 后 ， 按 照 进 化 状态 
情况 以 及 归档 集 规模 情况 ， 将 归档 集中 排序 靠 后 的 多 余 粒 子 裁 
BY fa Bl Hy 
3.2 全 局 最 优 粒子 的 选取 

粒子 群 全 局 最 优 粒子 决定 着 种 群 进化 方向 ， 也 影响 着 种 群 
的 选择 压力 和 种 群 的 多 样 性 。 因 此 ， 为 了 选取 综合 性 能 最 强 的 
粒子 作为 全 局 最 优 粒子 必须 建立 一 种 有 效 的 、 合 理 的 选取 策略 。 


at 


六 


本 文 利 用 旋转 基 可 视 化 技术 的 优点 ， 采 用 旋转 角度 和 旋转 基 折 
线 距 离 作 为 全 局 最 优 粒 子 的 选取 依据 。 多 越 小 ， 旋 转 基 折 线 距 
离 越 大 ， 该 粒子 的 多 目标 综合 性 能 越 好 ， 公 式 如 下 
3 cos Of 
“= arctan(4 ) 13 
> sin Of "9 


全 局 最 优 粒子 的 一 般 选 取 过 程 如 下 : 首先 按照 旋转 角度 排 
序 ， 选 择 旋转 角度 最 小 的 粒子 作为 全 局 最 优 粒 子 ， 如 果 有 若干 
粒子 的 旋转 角度 相同 ， 则 比较 其 旋转 基 折 线 距离 ， 然 后 选择 旋 
转 基 折 线 距 离 最 长 的 粒子 作为 全 局 最 优 粒子 。 

根据 以 上 分 析 可 知 ， 在 文献 [8] 算 法 的 基础 上 ， 对 归档 集 的 


粒子 fi f fo Lilols 格 占 优 强 度 旋转 基 角度 
Pı 0.6010 0.0003 0.0047 4 1 1 0 0.4636 
P 0.5270 0.0309 0.0232 4 4 1 0 1.0242 
Ps, 0.4715 0.0309 0.0655 3 4 2 0 1.0383 
Py 0.2081 0.0031 0.2002 1 1 4 0 0.5215 
其 次 ， 在 全 局 最 优 粒子 的 选取 中 ， 文 献 [8] 采 用 了 格 占 优 强 
度 和 密度 混合 算法 。 实 际 上 ， 格 占 优 强度 算法 在 一 些 特定 条 件 
下， 区 分 粒子 的 优 劣 能 力 不 足 。 例 如 ， 假 设 存在 4 个 粒子 3 个 
标 最 优 解 ， 其 格 占 优 强度 皆 为 0， 有 具体 如 表 1 所 示 。 为 此 ， 
必须 利用 密度 指标 进行 二 次 排序 才能 选 出 全 局 P1。 格 占 优 强度 


算法 的 分 辩 力 之 所 以 较 弱 ， 主 
旋转 基 的 可 视 化 选取 策略 ,可 以 直接 


要 是 


于 该 算 最 优 粒子 。 本 文 的 


区 分 4 个 粒子 的 优 劣 程度 ， 


序 计算 即 可 获知 最 优 粒子 为 法 在 算计 平行 格 坐标 时 采用 


KA 
向 上 取 整 的 方法 ,会 出 现 计算 误差 ,平面 网 格 划 分 的 格 数 越 少 ， 
损失 的 信 ， 


转 基 角 度 


旋转 基 角 度 时 没有 信息 损失 ， 


辨 能 力 ， 这 也 是 可 视 


4 


4.1 


Wr 


化 技术 的 


筷 越 多 ， 造 成 的 误差 也 就 越 大 。 相 反 ， 本 文 算法 的 旋 
决 射 函数 会 将 粒子 之 间 的 微小 差别 放大 。 另 外 ， 计 算 


因此 对 粒子 的 优 务 具有 较 强 的 分 


主要 特点 。 


rtMOPSO 算法 流程 与 时 间 复 杂 度 分 析 


算法 的 流程 


a) 种 群 参 数 初始 化 。 粒 子 群 种 群 主要 参数 包括 种 群 的 归档 
集 规模 天 、 种 群 的 规模 n、 最 大 循环 次 数 tmas AH 


Be A) 


子 


cl、cz。 将 归档 集 粒子 数量 k, EE ee t 等 参数 清 零 。 


b) 进 化 ; 


te 


AXPE; 然后 利 


后 ， 


按照 2.2 节 的 进化 状态 的 


天 态 判 断 。 将 归档 集 粒子 按照 式 (3) (4) 
J (7) (8) HES 


Ba 


归档 集 管理 与 


与 维 


子 。 


d) 根 据 3.2 节 策 略 更 新 全 


9) 按 照 式 (2) 更 新 粒子 群 


局 最 优 粒子 。 
种 群 。 


映射 到 旋 


RAS IB A FEHB; Bic 
值 判断 进化 状态 。 
佳 护 。 根 据 进化 状态 , 按照 3.1 节 的 规则 对 
归档 集 粒子 排序 并 剪裁 多 余 粒 


如 果 大 于 或 等 于 最 大 循环 次 数 则 停止 运行 ， 输 出 最 优 解 


R: 否则 ， 跳 到 b)。 
42 算法 时 间 复 杂 度 分 析 

在 rtMOPSO 算法 中 , 为 了 便于 计算 , 设置 粒子 群 种 群 和 归 
档 集 的 大 小 都 为 N， 多 目标 优化 的 目标 数量 为 MY， 粒子 维 数 为 
D. 根据 算法 的 流程 ，rtMOPSO 算法 主要 步骤 的 时 间 复 杂 度 的 
分 析 如 下 : 

Step1， 初 始 化 复杂 度 为 O(ND); Step2, 旋转 基 启 区 映射 为 


Chir 
录用 入 RAE, F: 基于 旋转 基 技 术 


线性 操作 , RARA Ze J EA INT Ti] 32 28 BE ONM); 在 Step3 中 ， 复杂 度 为 O(NlogN); Step5, 粒子 种 群 的 位 置 和 速度 更 新 为 线性 

首先 从 粒子 群 种 群 中 选择 有 效 粒 子 进入 归档 集 ， 其 复杂 度 为 ” 操作 。 综 上 所 述 ， 本 文 算法 的 最 大 时 间 复 杂 度 的 值 为 OMN?)。 

O(MN?)。 其 次 ， 执 行 维护 策略 时 需要 分 别 计算 收敛 程度 、 分 布 。 因此 ,本文 算 法 和 文献 [8] 算 法 以 及 主要 采用 Pareto 支配 .NSGA 

程度 和 辅助 判断 指标 , 其 时 间 复 杂 度 为 ONM), 多 指标 排序 的 ”快速 排序 等 方法 的 多 目标 粒子 群 优化 算法 的 时 间 复 杂 度 相同 。 
表 2 测试 结果 


Problem NSGAII SPEA2 MOEAD dMOPSO SMPSO tMOPSO 

DTLZ1 8.0280e-3 6.8079e-2 5.8709e-3 1.6370e+1 8.5590e+0 5.7273e-3 

DTLZ2 1.1326e-3 1.0332e-3 S.0483e-4 1.5412e-2 3.8664e-3 6.1321e-4 

DTLZ6 4.9899e-6 4.7047e-6 3.6114e-3 1.9039e-2 4.8176e-6 3.7941e-6 

on ZDT1 3.049 1e-5 2.3469e-5 8.6386¢e-5 6.0845e-4 3.8624e-5 0.1350e-5 

ZDT2 1.8500e-5 4,6469e-6 1.8572e-4 5.5901e-4 2.9975e-5 2.7368e-6 

ZDT3 3.9224e-5 3.2040e-5 2.1738e-4 4.1661e-4 4.3316e-5 2.5246e-5 

DTLZ1 4.855le-1 2.7821e-1 1.6938e-1 7.9282e-1 8.3284e-1 1.7325e-1 

DTLZ2 5.3215e-1 9.1660e-2 1.7041e-1 2.3881e-1 4.0660e-1 8.4357e-2 

aa DTLZ6 7.3369e-1 1.2879e-1 2.0492e+0 1.9683e+0 4,4293e-1 1.1258e-1 

oe ZDT! 4.4291e-1 1.2359e-1 2.7783e-1 2.7843e-1 4.1422e-1 1.0139e-1 

ZDT2 4.4236e-1 1.2982e-1 1.4822e-1 3.2305e-1 4.1902e-1 1.0846e-1 

ZDT3 4.5690e-1 1.1906e-1 5.8705e-1 5.6158e-1 3.9600e-1 0.9937e-1 

DTLZ1 1.2142e-1 1.3848e-1 1.2491e-1 5.3803e-3 1.3352e-2 1.8945e-1 

DTLZ2 7.0738e-1 7.3819e-1 7.4475e-1 5.5509e-1 6.7827e-1 7.0857e-1 

A DTLZ6 1.3270e-1 1.3316e-1 1.2103e-1 1.2121e-1 1.3284e-1 1.4183e-1 

pam ZDTI1 8.709 1e-1 8.7186¢-1 8.7036e-1 8.6267e-1 8.7034e-1 8.9032e-1 

ZDT2 5.3771e-1 5.3858e-1 5.3734e-1 4.4506e-1 5.3721e-1 5.5126¢e-1 

ZDT3 1.0202e+0 1.0237e+0 1.0139e+0 1.0101e+0 1.0232e+0 1.0328e+0 

表 3 部 分 对 比 算法 的 仿真 图 形 
函数 dMOPSO SMPSO NSGAII rtMOPSO 

ZDT1 
ZDT2 


ZDT3 
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DTLZ1 


DTLZ2 


DTLZ6 


5 ”实验 与 分 析 


51 测试 算法 及 对 比 算法 选择 
为 了 检验 改进 算法 rtMOPSO 的 算法 性 能 , 本 文选 取 ZDT05] 


系列 和 DTLZ09 系 列 中 6 个 函数 作为 测试 集 , 包 括 ZDT1、ZDT2、 


ZDT3 和 DTLZ1、DTLZ2、DTLZ6。 其 中 ZDT1、ZDT2、ZDT3 
三 个 函数 测试 改进 算法 在 凸 和 非 凸 以 及 不 连续 的 二 维 Pareto 前 
沿 上 寻 优 的 能 力 ; DTLZ1, DTLZ2, DTLZ6 三 个 函数 测试 对 多 
目标 、 多 变量 函数 的 寻 优 能 力 。 对 比 算法 选择 了 dMOPSONT、 


SMPSOU8 共 两 种 多 目标 粒子 群 算 法 以 及 NSGA-III9、SPEA2PO、 


MOEA/DP0 共 三 种 典型 的 进化 算法 进行 对 比 测试 。 
在 比较 仿真 实验 中 ， 为 保证 实验 的 客观 和 公平 性 ， 所 有 的 
测试 函数 的 主要 参数 设置 一 致 。 例 如 ， 粒 子 种 群 规模 为 100, 
归档 集 规模 为 100， 最 大 评估 次 数 为 30000; ZDT 系列 测试 的 
标 维 数 设置 为 2，DTLZ 系列 测试 的 目标 维 数 设置 为 3; 进化 
算法 的 交叉 率 设 为 0.9， 变 异 概率 为 0.1。 
5.2 评价 方法 
为 了 评价 对 比 算法 的 收敛 性 和 多 样 性 选择 以 下 三 种 指标 : 
a) 世 代 距 离 (generational distance, GD )22 反 映 获 取 的 Pareto 
解 集 靠近 真实 值 的 程度 。 其 值 越 小 ， 说 明 越 接近 真实 的 前 沿 。 
b) 多 样 性 指标 (spread ) 091 主 要 表明 对 解 集 分 布 均匀 情况 。 
其 值 越 小 ， 说 明 解 集 的 分 布 性 越 好 。 
c) 超 体积 (hypervolume，HV)JP3 可 以 衡量 解 集 覆 盖 的 目标 空 
间 区 域 范 围 并 同时 度量 了 收敛 性 和 多 样 性 的 性 能 。HYV 的 值 越 
大 ， 说 明 解 集 的 综合 性 能 越 好 . 
5.3 ”实验 结果 的 比较 与 分 析 


= 


= 


2 和 3 所 示 。 


其 中 表 2 是 将 rtMOPSO 和 各 种 对 比 算法 运行 30 


次 后 获得 实验 
试 实验 中 ， 在 


结果 。 根 据 表 2 可 知 ，rtMOPSO 在 二 维 目 标的 测 
GD, Spread 指标 上 都 优 于 其 他 对 比 算 法 ， 因 此 


对 2 维 目标 问题 ， 本 文 算法 表现 出 较 好 的 收敛 性 和 多 样 性 ， 对 


3 维 目标 问题 


DTLZ6 测试 ， 
对 于 大 部 分 算 


只 有 MOEA 的 DTLZ1、DTLZ2 的 性 能 和 本 文 算法 接近 ， 


， 总 体 上 来 看 ， 本 文 算法 的 收敛 性 也 比较 突出 ， 


本 文 算法 性 能 与 SPEA2 持平 。 实 验 数据 表明 ， 


法 而 言 , 对 DTLZ1 中 的 测试 结果 略 有 不 足 , 但 本 


文 算法 在 2、3 维 目标 的 测试 过 程 没 有 明显 缺陷 。 


为 了 进 
引进 行 比 对 


生 


步 说 明 本 文 算法 的 有 效 性 ， 可 通过 算法 的 对 综合 
。 本 文采 用 了 HV 指标 ， 同 时 从 收敛 性 和 多 样 性 
。 从 表 2 可知， 虽然 Spea2 和 MOEAD 是 比较 优 


= 


Ay Abe Se 


fi 
行 对 比 说 明 


有 较 好 的 综合 性 能 ， 但 本 文 算法 rtMOPSO 在 HV 


LIZ » 


进 
秀 
方 


表现 更 佳 
前 沿 ， 获 得 最 
为 了 更 加 


， 这 表明 本 文 提 出 的 算法 能 够 更 好 的 覆盖 Pareto 
优 解 集 。 
直观 、 清 晰 地 显示 本 文 算法 rtMOPSO 的 优越 性 ， 


在 表 3 中 ,本 文 给 出 了 rtMOPSO 和 部 分 对 比 算法 的 仿真 曲线 。 


根据 表 3 可 知 
足 之 处 。 例 如 
的 性 能 稍 显 不 


， 其 他 算法 在 三 个 性 能 指标 上 存在 不 同 程度 的 不 
dMOPSO、SMPSO 两 种 粒子 群 算法 在 3 维 目标 
E, 尤其 是 在 DTLZ1 的 测试 实验 中 , 解 集 的 分 布 


性 较 差 ， NSGAI 算 法 对 ZDT3 和 DTLZ2 的 测试 中 各 种 性 能 较 


差 。 本 文 算法 的 解 集 无 论 在 2 维 的 ZDT 真实 前 沿 曲线 上 ,还 是 


3 维 的 DTLZ 
， 本 文 算法 
综 上 所 述 


地 


真实 前 沿 曲 面 上 都 非常 接近 真实 解 。 因 此 综合 考 
优 于 5 种 对 比 算法 。 
， 本 文 算法 rtMOPSO 具有 较 大 优势 主要 原因 是 


本 文 算法 采用 


为 了 便于 实验 结果 ， 本 文 将 测试 结果 列表 说 明 ， 有 具体 如 表 


T, tem SA 


了 旋转 基 可 视 化 技术 ， 有 利于 决策 者 选择 最 优 粒 
档 集 的 管理 质量 并 确保 选取 了 全 局 最 优 粒子 。 另 
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外 ， 


中 
种 
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Pareto 前 沿 的 分 布 和 进化 状态 ， 有 利于 决策 者 选择 合适 的 策略 


在 粒子 群 算法 中 ， 采 用 了 进化 状态 监测 策略 。 在 进化 过 程 
判断 进化 状态 后 ， 选 择 合理 的 归档 集 的 排序 和 搜索 策略 ， 使 
群 同 时 具有 较 好 的 收敛 性 和 多 样 性 。 


tF 


结束 语 


网 化 技术 ， 


本 文 改 进 了 旋转 基 可 视 可 以 清晰 呈现 归档 集 粒子 


当 集 和 全 局 最 优 粒子 。 另 外 ， 借 助 Pareto KiF Ae i HY ABE 
监控 种 群 进 化 的 情况 并 根据 总 群 进化 的 状态 自动 选择 不 同 的 
策 


， 解 决 了 粒子 群 进化 过 程 中 快速 收敛 和 多 样 性 的 


引 


文 利用 旋转 基 可 视 化 技术 设计 了 双 排 序 的 档 
进 了 传统 的 按 密度 排序 的 单一 排序 方法 ， 使 
多 样 性 两 个 方面 同时 满足 要 求 。 最 后 ， 使 用 
转 基 角 度 和 距离 ， 提 出 了 全 局 最 优 粒 子 的 选择 策略 。 该 策略 
合 旋转 基 可 视 化 技术 的 优点 ,便于 选择 综合 性 能 最 强 的 粒子 ， 
导 种 群 向 Pareto 前 沿 靠近 。 
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